Tunable coupling between three qubits as a building block for a superconducting 

quantum computer 
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Large scale quantum computers will consist of many interacting qubits. In this paper we expand 
the two flux qubit coupling scheme first devised in [Phys. Rev. B 70, 140501 (2004)] and realized 
in [Science 314, 1427 (2006)] to a three-qubit, two-coupler scenario. We study L-shaped and line- 
shaped coupler geometries, and show how the interaction strength between qubits changes in terms 
of the couplers' dimensions. We explore two cases: the "on-state" where the interaction energy 
between two nearest-neighbor qubits is high, and the "off-state" where it is turned off. In both 
situations we study the undesirable crosstalk with the third qubit. Finally, we use the GRAPE 
algorithm to find efficient pulse sequences for two-qubit gates subject to our calculated physical 
constraints on the coupling strength. 



I. INTRODUCTION 

Quantum computers hold great promise to dramati- 
cally alter the way we perform computation. They can 
be constructed from qubits, two-level quantum systems 
which act as their fundamental building blocks. Al- 
though there are different approaches in actually building 
two-level quantum systems, in recent years superconduct- 
ing devices have seen some considerable success. There 
have been multiple proposals for designs of superconduct- 
ing qubits and of their possible coupling mechanisms 
\\^ . In this paper, we look at inductively coupled flux 
qubits along with dc Superconducting QUantum Inter- 
ference Devices (DC-SQUIDs) used for tunability. The 
DC-SQUIDs are treated as high excitation energy quan- 
tum objects and are assumed to always stay in their 
ground states [ll|, El- This method was first theoreti- 
cally described in Il3i . and subsequently demonstrated 
experimentally i n |l 4| . Similar coupling approaches were 
presented in 

This paper is organized as follows. Section HIl describes 
the system under study, in particular the different geome- 
tries that are explored. Section IIIII looks at the mathe- 
matics of coupling and describes how it varies in terms of 
physical system parameters. Next, Section [IVI shows nu- 
merical results when specific experimentally achievable 
parameters are used. Finally Section |V] looks at the ap- 
plication of the GRAPE algorithm when performing two- 
qubit gates. The conclusion is presented in Section [VIl 



II. SETUP 

A flux qubit can be constructed from a superconduct- 
ing loop interrupted by three Josephson junctions [l|. 
A current flowing in one direction can represent a com- 
putational basis quantum state |0) and in the opposite 
direction, a state The strength of the coupling be- 
tween qubits is mediated by a circulating current in a 




FIG. 1. Two different arrangements of flux qubits 1, 2, 3 and 
couplers A, B. "a" represents the arm length of the couplers, 
and "to" their width. We assume that the wires of couplers 
A and B overlap at most in two places around qubit 2. a) 
Straight-line shape, b) Capital "L" shape. 



DC-SQUID which is inductively coupled to each of the 
qubits. This circulating current can be controlled by the 
applied flux as well as the bias current, both of which can 
be tuned experimentally. Furthermore, the DC-SQUIDs 
may also be used for qubit readout [2(]| - |22| . 

The performance of any inductively coupled qubit ar- 
chitecture will strongly depend on the geometry of the 
qubit-coupler layout. The couplers themselves are based 
on the long-distance design proposed in [2^. In this pa- 
per we study two geometries, which can potentially serve 
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as building blocks for topologically motivated approaches 
to quantum computation [23, with high error thresh- 
olds [2B-I2II ■ In these approaches the physical qubits are 
typically laid out in a (for example) 2D square lattice. 
The geometries we study, which may be used to build 
such a lattice, are both shown in Fig. [1] In a), the cou- 
pling DC-SQUIDs form a straight line and in b), they 
enclose a 90° angle reminiscent of a capital letter "L". 
In both cases we have two qubits at the end points (la- 
beled "1" and "3") and the third (labeled "2") in the 
middle. The middle qubit is therefore surrounded by 
both of the couplers. We consider the qubits at the 
edge to be nearest- neighbours to the middle qubit. The 
qubits are placed inside the couplers, hence maximizing 
the qubit-coupler inductance. We take the edge length 
of the square qubits to be 44/im, roughly the same as 



those used in [IJ]. The distance between the qubit and 



the surrounding wire of the coupler is between 2 and 
4//m. To understand what impact qubit separation has 
on the interaction energy, we vary both the couplers' 
arm length and arm width (shown as "a" and "w" re- 
spectively in Fig. [T]). The arm length a varies between 
25/im and 300/im and wc study cases of the arm width 
w = 48^m and w = 2/im (although in the latter case 
we keep a 10 x IO/ito^ loop in the center of the arm to 
simplify delivering flux through the couplers. We stress 
that even though the coupling DC-SQUIDs are the same 
in both geometries, we expect the coupling energies to 
differ as even for the same length couplers, the distance 
between qubit 1 and 3 will be smaller in the 'L' shaped 
scenario. In order to numerically calculate the induc- 
tance for any given configuration, weperform finite el- 
ement simulations using FastHenry [29[ and assume we 
have aluminum wires with a penetration depth of 51nm 
[30| . The result is a 5 x 5 inductance matrix (see section 
IIVI for an example) . Since qubit 2 is always surrounded 
by two DC-SQUIDs, an experimental realization of this 
scenario would require crossing the wires of the couplers 
in at least two places. For simplicity we assume a fully 
symmetric situation meaning the qubits that are situated 
at the edges look the same to the qubit in the middle. 




FIG. 2. The three flux qubits labeled 1, 2 and 3, and two 
DC-SQUIDs labeled A and B, form a three-qubit, two-coupler 
system. 7„,i is the phase diff'erence across the ith junction in 
a coupler n. and J„ are the applied flux and circulating 
current of coupler n respectively. lo is the critical current 
of the Josephson junctions used in the coupling DC-SQUIDs 
— for simplicity we assume that all the junctions, and hence 
their critical currents, are the same. 



III. THEORY 

In this section we expand on the results presented in 
[isj . First, we define the parameters of our system of 
which a schematic representation is shown in Fig. [51 Each 
variable j„^i is the phase difference across the ith junc- 
tion in coupler n. $„ is the applied flux threaded through 
coupler n, /„ its bias current and J„ its circulating cur- 
rent, /q is the critical current of the Josephson junctions 
used in the coupling DC-SQUIDs — for simplicity we 
assume that all the junctions, and hence their critical 
currents, are the same. We further assume that we are 
able to control the flux independently through each of 
the DC-SQUIDs and qubits. This may need multiple 
bias coils and careful calibration to cancel cross-control 
effects. We can write an effective Hamiltonian which de- 
scribes the jth flux qubit [l| as 



(3.1) 



where (Tz and Ux are Pauli matrices, ej is the energy bias 
which depends on the flux threaded through the qubit 
and Aj is the tunneling energy, which we take to be fixed 
at fabrication time and usually of the order of a few GHz. 
We stress that modified flux qubit designs, with variable 
Aj have been proposed [l[ and demonstrated experimen- 
tally [llHsJl, but in principle one only needs control of 
one of the qubit's quadratures to do arbitrary single- 
qubit gates. Next, we note that each qubit i interacts 
with any other qubit j via the direct mutual inductance 
as well as via the mutual inductance mediated by the 
couplers. We can therefore write the total Hamiltonian 
of the system as 



H 



j = 1.2,3 



(3.2) 



where Hj are the Hamiltonians of the individual qubits, 
and Kij correspond to the coupling energy between 
qubits i and j 
as 



Each Kij can be further described 



E 

n=A,B k=A,B 



dJn 



(3.3) 



with M the inductance matrix of the system. We eas- 
ily see that the first term corresponds to the direct cou- 
pling between the qubits, and the other terms are in- 
volved in the indirect interaction through the coupling 
DC-SQUIDs. The negative sign in front of the second 
term is justified in [3a, H^- Each of the couplers n (for 
n = A, B) can be controlled by either applying a bias 
current /„ and/or a bias fiux By adjusting these 

two parameters we will study how the coupling between 
any pair of qubits can be either enhanced or reduced. 
The interaction energy in Eq. 13.31 involves the transfer 
functions of the form ^ which wc have to calculate. 

We start with Kirchhoff's current laws for coupler n. 
In the slowly varying limit of the junctions (l3| . we can 



3 



write 



/„ = lo sin 7„i + lo sin 7„2 
2 J„ = /o sin 7„i - Jo sin 7„2 ■ 



Next, using jn- — ''"^ ^''"^ and -f„ 



2/osin7„+ cos7„ 



7nl+7Ti2 



rewrite Eqs. 13.41 as 
and 

J n 



/osin7„_ cos7„+. 



(3.4) 



(3.5) 



(3.6) 



From flux quantization [371 138| , we have another con- 
dition, which relates the flux passing through the DC- 
SQUID to phase differences across the junctions. This 
lets us write 



Jn as well as the phases 7„_ and 7n+. Next we can 
implicitly differentiate Eqs. 13.51 13.61 and 13.71 in order to 
obtain the transfer functions and (here k is just 
a coupler index with k ^ n), which we will need to cal- 
culate the interaction energy Kij. Noting that the bias 
current does not depend on the flux, and differentiating 
Eq. 13.51 gives 



= (sin7„+cos7„_) (3.8) 



9$ 



which further leads to 



dln+ .2 .2 9ln- 

— — = tan 7„+ tan 7„_ — — 



(3.9) 



In 



MnkJk)- 



(3.7) 



By fixing the applied flux and bias current /„ we can 
solve Eqs. 13.51 13.61 and 13.71 to obtain circulating currents 

I 



In total, this gives us four equations, one for every com- 
bination oi n,k = A, B. Differentiating Eq. 13.61 and 
substituting in Eq. 13.91 leads to 



dJn 



Iq (cos7„_ cos7„+ - sin7„_ sin7„+ tan7„_ tan7„+) 



(3.10) 



Finally we are left with obtaining which will de- Substituting the results of Eqs. 13.121 and 13.111 into 13.101 

pend on whether n and k are the same or different. Hence leads to 
differentiating Eq. 13.71 and taking n = k gives 



djn- TT / , , dJn . , dJk 



and when n ^ k 



TT / dJn dJk 



9$, 



(3.11) 



(3.12) 



1 - tan^ 7„+ tan^ 7„_ 



\ + tan^ 7„+ tan^ 7„_] 



1 - A/„ 



dJk 



(3.13) 



tan^ 7fc+ tan^ ^k- 



dJk ^ ^ 

d-^n 2Li 1 + ili^ [1 - tan2 7^+ tan^ 7^.] 



dJn 



(3.14) 



r 



where we have taken 



$0 



-n „ . (3-15) 

27r7o cos7„_ cos7„+ 

to be the Josephson inductance. We easily observe that 
in the case of one coupler, M„ ^ = in Eq. 13.131 and Eq. 
13.141 does not come into play. This reduces the result to 



the one obtained in the original paper [l3|. Hence, for a 
given geometry, applied fluxes <^a,^b and bias currents 
Ia,Ib, we now have a means of calculating the Hamilto- 
nian from Eq. 13.21 
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IV. COUPLING RESULTS 



In our numerical calculations, we take the qubits' per- 
sistent current values to be Ip = 0.46/iA and the critical 
currents of the Josephson junctions that are used in DC- 
SQUIDs as 0.11/iA — both in the experimentally achiev- 
able ranges. In order to answer whether all couplings can 
be turned off, or whether we can selectively make certain 
interactions strong, while keeping others weak, we need 



to understand what happens to all the coupling terms 
for all possible input parameters. To do this, we scan 
both of the applied fluxes and <f>B between — (f>o/2 
and $0/2 and the bias currents between and values 
approaching the critical currents of the couplers (which 
depend on the threaded fluxes). The inductance matrix 
is calculated for each geometry and used to obtain values 
for all possible Kij. For example, in the L-shapcd geom- 
etry with the thin arm width oi w — 2/im and the arm 
length of a = 50/xm, its value is 
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where the columns (and rows) correspond to the first, 
second, third qubits, coupler A, and coupler B respec- 
tively. 

Fig. [3] shows the coupling strength of the L-shape ge- 
ometry, where the arm width "w" is 2/xm and the arm 
length "a" is 100/im. The horizontal axis corresponds to 
the values of coupling energies if 12, the vertical to if23 
and finally the "out-of-page" direction to the crosstalk 
term Ki^. The plot was obtained by varying the fluxes 
and biasing currents through both couplers. Thousands 
of data points were calculated, which were then used to 
obtain a surface that spans different combinations of the 
three coupling energies. The process was then repeated 
for varying arm lengths and arm widths of the coupling 
DC-SQUIDs. From the results one can note that the 
crosstalk term A'13 can be treated as a single- valued func- 
tion of A'i2 and ii'23- We will use this fact in the next 
section. 

The first question we wish to address is whether we can 
keep the coupling between two nearest-neighbour qubits 
(say 1 and 2) high, the other nearest-neighbour coupling 
(say 2 and 3) turned off, while eliminating the crosstalk 
(interaction between qubits 1 and 3). We concentrate on 
this scenario as we envision doing gates between nearest- 
neighbors (here qubits 1 and 2), and would like the qubits 
that are not involved in the gate (qubit 3) to be disturbed 
as little as possible. In practice, we will find that it is 
not possible to completely turn off the interactions with 
the third qubit, and hence all the values of K will need 
to be varied, but during a gate between qubits 1 and 2 
the interaction 7^12 will dominate. We can get the an- 
swer to this first question by studying plots like the one 
presented in Fig. [3] We first fix K23, (vertical axis) at 
zero, and then traverse in the horizontal direction along 
the dashed line. Figs. |3] and [S] show cross-sections for 
the line-shaped and L-shaped geometries, with the arm 




O.Q 0.5 LO 

Kn (GHz/h) 



FIG. 3. (Color online). Interaction energies 7^12 (horizontal 
axis), K23. (vertical axis) and Js'ia (out-of-page direction) for 
an L-shape geometry with couplers' arm width of ui = 2^m 
and arm length a — lOO^m. In the calculations we take /i = 
0.46M and Jo = O.IIaiA. 



width w ~ 2fim. In the line-shaped geometry scenario 
(Fig. 2]), with shortest considered arm length a = 25/im 
and K23 tuned to approximately 0, we find the crosstalk 
Ki3 to be 47MIIz. The crosstalk only varies slightly as 
we increase the coupling Ki2. In the L-shape geometry 
(Fig. [5]), the situation is less favorable, as expected, since 
qubits 1 and 3 arc closer together and their direct mutual 
inductance is larger. Furthermore, they interact with the 
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coupler that is further away from them more than is the 
case in the line-shaped geometry. From Fig. [5] we see that 
the best we can do given short couplers with a = 25/im 
is i^ia 160MHz. For concreteness, if one considers a 
"high" value of coupling 1-2 to be « l.OGHz, these re- 
sults lead to the unwanted crosstalk of 5-16% percent 
when the qubits are closest together. However, as ex- 
pected and shown in Figs. [5] and SI we can reduce the 
unwanted crosstalk term iiTis by simply increasing the 
arm length of the coupling DC-SQUIDs — for example 
in a case of a = 300/im, for A'12 > l.OGHz, the magnitude 
of Ki^ is smaller than 2MHz in all the geometries. We do 
stress however that one cannot keep increasing the arm 
length indefinitely, as the process also decreases the max- 
imum achievable energy that can be mediated between 
nearest-neighbours qubits. We have further found that 
increasing the arm width of the couplers to w = 48/ini 
made no qualitative difference in the results above, how- 
ever for the same set of parameters the magnitudes of the 
interaction energies Kij were smaller. This is due to the 
fact that the mutual inductance elements between these 
wider coupling DC-SQUIDs and qubits are smaller than 
in cases of the narrow couplers. The discussion above 
leads us to conclude that minimizing the crosstalk cou- 
pling (term /V13) can be done to a great degree, even 
when fairly small DC-SQUIDs are used and more impor- 
tantly in both of the studied geometries. 
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FIG. 5. (Color online). The size of the crosstalk interaction 
K\i as a function of K12 in the L-shaped geometry scenario 
where the coupling K2i has been taken to be zero. In the 
calculations we take Ii = I2 = Is ~ and Jo — 0.11/.tA. 



the results in Fig. [6] for both geometries and two dif- 
ferent arm widths. From the data it is clear that there 
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FIG. 4. (Color online). The size of the crosstalk interaction 
A"i3 as a function of K12 in the line-shaped geometry scenario 
where the coupling K23 has been taken to be zero. In the 
calculations we take h = I2 = I3 = 0.46/iA and lo = 0.11/iA. 

The next question we wish to answer is whether the 
interaction between all pairs of qubits can be turned off 
completely. This may be useful when one wants to per- 
form single gate operations on any (or all) of the qubits 
without affecting the others. To determine this, we once 
again turn to Fig. [3] (and its analogues for other geome- 
tries and values of arm lengths and widths — not ex- 
plicitly shown). We can now look at the crossing of the 
two dashed lines, which corresponds to K12 = K23 = 0, 
and figure out the value of the crosstalk (ATia). We plot 



-0.1 - 



-0.15 



-0.2L 



□ Line shape, w=2[im 
-e— Line stiape, w=48iim 
■ ► - L shape, w=2\im 
— L shape, w=48^im 



50 



100 150 200 
Arm length (|im) 



250 



300 



FIG. 6. (Color online). A plot consisting of values for the 
crosstalk K13, while keeping K12 = K23 ~ 0. The horizon- 
tal axis represents different arm lengths of the coupling DC- 
SQUIDs. We find that the magnitude of K13 is substantially 
smaller in the case of the straight-line geometry, but even in 
the L-shaped scenario dies olT quickly. In the calculations we 
take 7i = /a = 73 = 0.46^, Jo = O.llAiA. 

is no complete off-state. The straight-line geometry is 
much less (over 3 times in a case of shortest arm length 
we studied) susceptible to the unwanted crosstalk if 13. 
However the crosstalk energy dies off quickly, and in a 
case of arm length a > 200/im, its magnitude stays below 
5MHz. Given that the typical energies of the qubits are 
of the order of a few GHz, a crosstalk strength of roughly 
three orders of magnitude smaller may not be very dam- 
aging. We also find that the width of the couplers has 
little effect on the on-off state. We further note that the 
crosstalk and off-state fidelity could in principle be im- 
proved by additional superconducting screening elements 
[391, which are not examined in our work. Others have 
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also shown that there is a distinct crosstalk mechanism 
related to the breakdown of the adiabatic approximation, 
discussed, e.g., in [l^l and exemplified in [40|. We chose, 
however, in this work to explore control that tolerates 
the crosstalk at hand. 



V. LOGICAL OPERATIONS 

A crucial component of any quantum computer is its 
ability to perform logical operations. These operations 
can be constructed from only single-qubit rotations and 
(at least) one two-qubit gate (such as a CNOT) [4l| . In 
[l3| it was shown how a two-qubit gate can be performed 
using results from (4^ . Here, in the three-qubit case, we 
concentrate on a different approach. We use numerical 
optimization via the GRAPE algorithm [i^ to find a 
collection of pulses that maximize the gate fidelity Fg, 
with 

Fg = Tr{[/t[/^.}2 (5.1) 

where U is the desired gate, and Un a calculated numer- 
ical approximation. GRAPE works by dividing the gate 
time into time-slices (often called time pixels) of constant 
pulse amplitude. The amplitudes at any time pixel are 
found by calculating the gradient of the fidelity function 
with respect to some well defined controls. Since the gra- 
dient points in the direction of increasing Fg^ over many 
iterations, one can converge to a pulse sequence that is 
(arbitrarily) close to the intended gate. We look at two 
gates, the first, a "natural" candidate for this type of 
system is the ViSWAP, and the other a CNOT gate. We 
stress that a CNOT gate is often used to describe the ac- 
tive error correction procedures in the topologically mo- 
tivated schemes from [H, [H, H!] . 




time (ns) 



FIG. 7. (Color online). A pulse sequence for a \/iSWAP 
gate between qubits 1 and 2. Since the system coupling is 
"natural" for this gate, no single-qubit rotations are needed. 
The pulses lead to the desired gate with 99.9% fidelity. In the 
calculations we take Ai/h = 5.0 GHz, A2/h = 5.4 GHz and 
As/Zi = 5.8 GHz. 



We start with the system Hamiltonian shown in 
Eqs. 13.11 and 13.21 and note that each ej is time depen- 
dent and can be written as 

ej{t) = 4e^(i) cos(wjt) + 4eJ(t) sm{ujjt). (5.2) 

Having both ej(t) and Cj{t) present, will later give us 
the ax and ay quadratures. We use both of these when 
creating CNOT pulses. We can now switch basis, to the 
eigenbasis of ax, namely {|+), |— )} and furthermore go 
into a rotating frame defined hy U = Ui ^ U2 ^ U3 with 
Uj — exp{—i^ait). We then assume that the qubits 
can be driven at their respective resonances (meaning 
Aj = Uj) and neglect the fast oscillating terms (do a 
RWA). AH this leads to an effective Hamiltonian in the 
rotating frame and {|+), |— )} basis 

He=Y.e^{t)ai-e^^{t)al 

i 

- Ki2{t)e-'^^'-^^^' (|100)(010| + |101)(011|) + h.c. 

- K23it)e-'^^^-^'^' (|001)(010| + |110)(101|) + h.c. 

- Ki3{t)e-'''^'-^'^* (|100)(001| + |110)(011|) + h.c. 

(5.3) 

We explicitly make the Kij terms time dependent to 
stress that they will vary. For demonstration purposes, 
in the case of yiSWAP, we chose the qubit energies to 
be Ai/h = 5.0 GHz, A2//1 = 5.4 GHz and Ag/Zi = 
5.8 GHz, and when creating the pulses for the CNOT 
gate we take all qubits to be on resonance — meaning 
Ai/h = A2/h = Aj,/h = 5.0 GHz. 

We constrain GRAPE by only allowing values of Kij 
which can be achieved using physical controls (as de- 
scribed in section ITV)). We can think of this as staying on 
a surface spanned by the triplet {K12, K2Z, K13). From 
Section ITVl we note that K13 can be treated as a single- 
valued function of Ki2 and K2Z. This slightly simplifies 
our procedure of finding the gradient of Fg, since we are 
only dealing with two independent variables that account 
for all the interactions in the problem. We further allow 
GRAPE to vary single-qubit controls (which effectively 
means ej(t) and ej(t) from Eg. 15. 3p . 

The resulting pulses for the %/ iS WAP and CNOT gates 
are presented in Figs. [7] and [5] respectively. Only the non- 
zero controls arc explicitly shown. We study the case 
of the L-shaped geometry with the arm length of 25/^m 
where the unwanted crosstalk coupling K13 is the largest 
(and as can be seen from the plots, never zero). 

In the case of viSWAP, it is clear that only coupling 
1-2 has to be varied, and no single-qubit rotations are 
needed on any of the qubits. Since we have chosen the 
qubits to be off-resonance with each other, in the rotating 
frame the terms K23 as well as the crosstalk K13 oscillate 
at frequencies (A2 — A3) and (Ai — A3) respectively. 
Because they are not too large and average out to zero 
over the length of the pulse, little correction is needed in 
the generated controls. 
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The situation is more complicated in the case of the 
CNOT gate, partly because it is a "non-natural" gate 
for the Hamiltonian we have, but also due to the fact 
that here we chose to keep all qubits on-resonance with 
one another. The resulting pulses require single-qubit 
controls and in particular both Ox and Oy quadratures on 
qubit 3 need to be varied. In contrast to the off-rcsonant 
case, to achieve high fidelities, qubit 3 must be repeatedly 
flipped and rotated to counteract the crosstalk term. 

In both cases, the pulses found perform with a fi- 
delity of 99.9%. That number however can be often in- 
creased (in which case the calculation can take substan- 
tially longer). The pulses in general are not unique and 
can depend on the initial guesses which GRAPE then 
modifies. 
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FIG. 8. (Color online). A pulse sequence for a CNOT gate 
between qubits I and 2. Single-qubit rotations are applied to 
all three qubits. For timely convergence of the GRAPE algo- 
rithm, we modulate both cr^ and Oy quadratures on the third 
qubit. Furthermore, we allow all three interaction energies to 
vary. The pulses lead to the desired gate with 99.9% fidelity. 
In the calculations we take Ai//i = 1^2 /h — A^/h = 5.0 GHz. 



Finally, let us note that we have shown how one has 
to vary the energies in the system to obtain the desired 
pulses. Alternatively one could concentrate on work- 
ing with values of physical parameters (currents, fluxes) 
which might be more useful when trying to devise an 
experimental demonstration of this system. 



VI. CONCLUSION 

In this paper, we discussed a system which may pro- 
vide a candidate for a many-qubit processor built out 
of superconducting circuits. In particular we extended 
the coupling mechanism first presented in [isj to a three- 
qubit, two-coupler scenario. We studied multiple geome- 
tries of the coupling DC-SQUIDs; a line-shape, and an 
L-shape, with couplers' arm lengths between 25/^m and 
300//m, and arm widths of 2 and ASfim. We looked 
at how much unwanted crosstalk (Ki^) there is while 
the nearest-neighbour interaction {K12 or K23) is high 
(w IGHz). With our experimentally achievable physical 
parameters, when the arm lengths are shortest, we found 
the maximum crosstalks of 47MHz and IGOMHz for the 
line-shape and L-shape geometries respectively. We also 
characterized the crosstalk in the off-state and found that 
for arm lengths greater than 200/im it was smaller than 
5MHz, which is three orders of magnitude smaller than 
typical qubit energies. Finally, we showed how one can 
use optimal control methods (GRAPE) to perform logical 
operations in the presence of crosstalk (while accounting 
for the constraints on all the coupling energies), and in 
particular we showed potential pulse sequences for the 
ViSWAP and the CNOT gates. 
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